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Abstract 



Since electronic and magnetic properties of many transition-metal oxides can be efficiently controlled by external factors such as 
C^) the temperature, pressure, electric or magnetic field, they are regarded as promising materials for various applications. From the 
C^l viewpoint of electronic structure, these phenomena are frequently related to the behavior of a small group of states close to the 
C^i Fermi level. The basic idea of this project is to construct a low-energy model for the states near the Fermi level on the basis of first- 
D principles density functional theory, and to study this model by modern many-body techniques. After a brief review of the method, 
the abilities of this approach will be illustrated on a number of examples, including multiferroic manganites and spin-orbital-lattice 
■^j- coupled phenomena in (R being the three-valent element). 
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1. Basic Idea, Purpose, and Methods of Realistic Modeling 

The oxide materials can be rather complex. Nevertheless, in 
many cases their electronic and magnetic properties are con- 
trolled by a small group of states located near the Fermi level 
and well isolated from the rest of the spectrum. A typical 
example of the electronic structure of YVO3, obtained in the 
local-density approximation (LDA), is shown in Fig. [TJ In this 
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Figure 1: (Left panel) Electronic structure of orthorhombic YVO3 in LDA. 
J.^ The shaded area shows the contributions of the 3d states of V. (Right panel) 
d Enlarged behavior of f9 g bands computed from LMTO basis functions (solid 
curves) and downfolded bands (dot-dashed curves). The corresponding bands 
in the left panel are shown by arrows. The Fermi level is at zero energy. 

case, the active states are twelve t2 g bands of predominantly V- 
character, which are separated by finite energy windows from 
other states, both from below and from above. Such a division 
of the entire electronic structure into the "active" and "inactive" 
(or low- and high-energy) states opens a formal way for com- 
bining the first-principles calculations, based on the density- 
functional theory (DFT), with the many -body treatment of some 



effective model, formulated rigorously in the restricted Hilbert 
space of "active" states. This is the main idea of realistic mod- 
eling of complex oxide materials. The purpose of this project is 
twofold: 

(1) To incorporate the physics of Coulomb correlations, 
which is greatly oversimplified in conventional LDA; 

(2) To provide a transparent physical picture for electronic 
and magnetic properties of complex compounds. In this sense, 
the realistic modeling can be regarded as supplementary ap- 
proach to more conventional first-principles electronic structure 
calculations, which are currently on the rise. 

Thus, the first step of the project is construction of the low- 
energy model (typically, the multiorbital Hubbard model) for 
the states near the Fermi level: 
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which would include the effect of other ("inactive") states in the 
definition of the model parameters of the Hamiltonian ([TJ. All 
these parameters are derived totally from the "first-principles" 
on the basis of DFT. We would like to emphasize that, although 
the derivations are inevitably based on some approximations, 
we do not use any adjustable parameters apart from these ap- 
proximations. The procedure of constructing the model Hamil- 
tonian was described in details in the review article fljj] . Briefly, 
in order to derive the one-electron part {t RR ,), we use the gener- 
alized downfolding method. For the isolated low-energy bands, 
this procedure is exact, as it is clearly seen in Fig. [TJ from the 
comparison of the original band structure, obtained in the lin- 
ear muffin- tin orbital method (LMTO) [2], and the one after 
the downfolding. The parameters of screened Coulomb inter- 
actions {U a p y s) are typically obtained by combining the con- 
strained DFT technique [3] with the random-phase approxima- 
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tion [4]. The latter is very efficient for treating the screening of 
correlated electrons by themselves [5], which can dramatically 
reduce the effective Coulomb interactions in the low-energy 
bands. For example, the bare Coulomb repulsion between 3d 
electrons is about 20-25 eV. However, for the low-energy bands 
in solids this value is typically reduced till 2-4 eV QJJ] . 

Once the model is constructed, it can be solved by using var- 
ious many-body techniques. In the present work, we typically 
start with the mean-field Hartree-Fock (HF) approximation, and 
take into account the correlation interactions by considering the 
perturbation theory expansion near the HF ground state. This 
procedure may be justified, if the degeneracy of the ground state 
is lifted by the crystal distortions QJJ] . 

Below, we present examples of realistic modeling for two 
types of transition-metal oxides. 

2. Spin-Orbital-Lattice Coupling in vanadates KVO3 

The vanadates (where R is the three-valent, typically 

rare-earth, element) have attracted a considerable experimen- 
tal and theoretical attention. All these compounds crystallize 
in the distorted perovskite structure. For example, considered 
in the present work YbV03 can have orthorhombic and mono- 
clinic modifications, realized at T- 15 and 75 K, respectively. 
However, a relatively small change of the lattice parameters 
may cause a dramatic reconstruction of electronic, and asso- 
ciated to it, magnetic structure. For example, at T- 15 and 
75 K YbV03 forms the so-called G- and C-type antiferromag- 
netic (AFM) structure, respectively (Fig. [2]). The phenomenon 




Figure 2: Distribution of the charge densities associated with the occupied fio 
orbitals (the orbital ordering) realized in the orthorhombic (left) and monoclinic 
(right) phase of YbVO^ in the HF approximation. Different magnetic sublat- 
tices associated with two opposite directions of spins in the AFM structure are 
shown by different colors. 

called "spin-orbital-lattice coupling" and is typically regarded 
as the test for various theories of the electronic structure. First 
attempts of realistic modeling of RWO3 are summarized in [Q]]. 
After that we have performed a systematic study for the entire 
series RW0 3 (R= La, Ce, Pr, Nd, Sm, Gd, Tb, Ho, Yb, Lu, and 
Y), using the experimental crystal structure available in the lit- 
erature Ja01- For all considered compounds, the effective Hub- 
bard model was constructed for ti & bands (Fig.[TJ in the basis of 
three Wannier orbitals per each V-site. The details will be pub- 
lished elsewhere. Here we present an example of the "canon- 
ical behavior" realized in YbV03, where the crystal distortion 



quenches the orbital structure in some particular configuration 
(Fig.|2j, which uniquely defines the type of the magnetic ground 
state, basically via Goodenough-Kanamori rules. In this case, 
once the degeneracy of f2 ? -levels is lifted by the crystal distor- 
tion, the correct type of the magnetic ground state can be suc- 
cessfully reproduced already at the level of HF approximation. 
For example, the stabilization energy is clearly the largest for 
the experimentally observed G- and C-type AFM states when 
YbVOi crystallizes in the orthorhombic and monoclinic struc- 
ture, respectively (Fig. [3}. Moreover, the correct AFM ground 
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Figure 3: Stabilization energies of the main AFM states in YbVC>3 relative to 
the ferromagnetic state as obtained in the HF approximation (dark blue area) 
and after taking into account the correlation interactions in the second order of 
perturbation theory (light blue area) and in the t-matrix theory (hatched area). 

state is additionally stabilized by correlation interactions, which 
are taken into account via perturbation-theory expansion near 
the HF solutions. We have considered two such technique for 
the total energy [8]: one is the regular second-order perturba- 
tion theory and the other one is the t-matrix theory. Both of 
them provide very consistent explanation for YbVO^, although 
the correlation energies obtained in the t-matrix theory for the 
AFM states are systematically smaller due to the higher-order 
correlation effects, which are included to the t-matrix, but not 
to the second-order perturbation theory. 

3. Inversion Symmetry Breaking in Manganites /?Mn03 

The multiferroic manganites, such as BiMnOi and TbMn03, 
are currently under very intensive investigation. The multifer- 
roicity means that the magnetic order in certain system without 
the inversion symmetry coexists with some finite ferroelectric 
polarization. The coupling of these two order parameters pro- 
vides a unique opportunity to control the magnetic properties 
by applying the electric field and vice versa. From this point of 
view, the most interesting is the situation, where the inversion 
symmetry is broken by the magnetic degrees of freedom. In 
order to study the microscopic origin of the inversion symme- 
try breaking in manganites, we have constructed the effective 
model in the basis of three tig and two e g orbitals at each Mn- 
site iSlllOt]. The main results can be summarized as follows. 

The mechanism of the inversion symmetry breaking is re- 
lated to the behavior of interatomic magnetic interactions, 
which depends on the distribution of occupied e g orbitals 
(Fig. 3|. Although the orbital ordering differs dramatically 
in monoclinic and orthorhombic manganites, the basic idea is 
rather generic: the nearest-neighbor interactions Jnn always 



2 



BiMnO, 




LaMnO, 




Figure 4: Typical distribution of the charge densities associated with the oc- 
cupied e g orbitals (the orbital ordering) realized in the monoclinic (BiMn03, 
left) and orthorhombic (LaMn03 , right) manganites with the notations of main 
magnetic interactions. 



compete with some long-range AFM interactions Jir. The ex- 
istence of Jlr is related to the fact that, due to the screening, the 
on-site Coulomb repulsion U is not particularly large (about 2.2 
eV in manganites). Therefore, other mechanisms of magnetic 
interactions, besides conventional superexchange (being of the 
order of 1/U), can be also operative. These mechanisms, which 
are of the higher orders than 1/U, together with the form of 
the orbital ordering, naturally explain the large values of Jm in 
some particular bonds H. When Jm becomes larger than /^v. 
they can lead to the formation of complex magnetic structures 
with the broken inversion symmetry. In monoclinic BiMnC»3, 
this is the collinear AFM structure, where each of the two sub- 
lattices (1,2) and (3,4) are coupled antiferromagnetically (see 
Fig. [4] for the notations of Mn-sites) ifToh . In TbMn0 3 , the 
ground state is the incommensurate spin-spiral, where the direc- 
tions of spins in the ab-planes vary as eR=(cosqR, sinqR, 0), 
and the spin moments between the planes are coupled antifer- 
romagnetically. Such a magnetic structure can be easily calcu- 
lated by using the generalized Bloch theorem, which combines 
the lattice translations with the spin rotations fiUl . The total 
energy minimum corresponds to the propagation along the or- 
thorhombic b-axis, in agreement with the experiment. How- 
ever, the obtained value ^=0.675 (in units n/b, Fig- EJ, is sub- 
stantially larger than experimental q/,= 0.25 [ 12] and 0.28 H 1 311 . 
reported for magnetic structures with the moments lying in the 
ab- and bc-plane, respectively. 

In order to resolve this discrepancy, we consider the relativis- 
tic spin-orbit interaction (SOI). In BiMnC<3, the SOI is respon- 
sible for the spin canting away from the collinear AFM struc- 
ture and formation of the net ferromagnetic moment. Thus, the 
ferroelectric response in BiMn03, caused by the hidden AFM 
order, coexists with the ferromagnetism and can be controlled 
by the magnetic field ifioll . The magnetic interactions of the 
relativistic origin play an important role also in TbMn03 iTfl . 
In the search for the true magnetic ground state of TbMn03, 
we have investigated several magnetic structures with different 
periodicity. The structure corresponding to the lowest energy is 
shown in Fig. [6] The periodicity of this structure is described by 
qb-0.25, in agreement with the experiment [12]. Nevertheless, 
it is no longer the uniform spin-spiral. It would be interesting 
to check our finding experimentally. 
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Figure 5: Dependence of the total energy (measured in meV per one formula 
unit) on the homogeneous spin-spiral vector q obtained in the Hartree-Fock 
approximation for TbMnOj without the relativistic spin-orbit coupling. 




Figure 6: Distribution of magnetic moments in the ab-plane obtained in the HF 
approximation with SOI for TbMnC>3 . The angles formed by Mn-moments in 
the bonds Mn-O-Mn are specified by numbers. 



4. Conclusions 

The realistic modeling combines the accuracy and pre- 
dictable power of first-principles electronic structure calcula- 
tions with the flexibility and insights of the model analysis. The 
first applications of this method are very encouraging. We hope 
that these ideas will continue to develop to become a power- 
ful tool for theoretical analysis of complex oxide materials and 
other strongly correlated systems. 
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